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MOORE-GIBSON-THOMPSON EQUATION WITH 
MEMORY, PART I: EXPONENTIAL DECAY OF ENERGY 

IRENA LASIECKA AND XIAOJUN WANG* 


Abstract. We are interested in the Moore-Gibson-Thonipson(MGT) 
equation with memory 


TUttt + autt + Au + bAut — g{t — s)Aw(s)ds = 0. 

Jo 

We first classify the memory into three types. Then we study how a 
memory term creates damping mechanism and how the memory causes 
energy decay. 
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1. Introduction 

We are interested in the Moore-Gibson-Thompson(MGT) equation with 
memory 


( 1 ) 


TUttt + oiUtt + c^Au + bAut — / g{t — s)Aw{s)ds = 0 

Jo 


where r, a,c^,b are physical parameters and A is a positive self-adjoint oper¬ 
ator on a Hilbert space H. The convolution term g{t — s)Aw{s)ds reflects 
the memory effects of materials due to viscoelasticity. Here the convolution 
kernel g satisfies proper conditions exhibiting “memory character”, which 

1 
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will be explained later. The form of w classifies the memory into the follow¬ 
ing three types 


(2) 

w[s) = u(s), 

(Type 1) 

(3) 

rc(s) = ut{s), 

(Type 2) 

(4) 

w{s) = Au(s) -|- ut{s). 

(Type 3) 


This is the first part of an ongoing project, in which a third order (in 
time) equation with memory is investigated. More specifically, we study 
how (and whether) a memory term creates damping mechanism and how 
(and whether) it causes energy decay. 

1.1. Background. Moore-Gibson-Thompson (MGT) equation arises from 
modeling high frequency ultrasound waves. Without memory, the linearized 
MGT equation reads 

(5) TUttt + autt + c^Au + bAut = 0. 

Gertainly this equation is in abstract form and it has a simple prototype 
where A = —A. In m, the well-posedness of ([5|) and the uniform decay 
of its energy is studied under proper functional setting and initial-boundary 
conditions. We shall revisit the results of m in Section 2. Spectral analysis 
for this model has been carried out in [23], which confirms the validity and 
sharpness of the results in m- 

A linear MGT equation is the prelude to nonlinear ones. The classical 
nonlinear acoustics models include the Kuznetsov equation, the Westervelt 
equation and the KZK equation, etc. This research field is highly active due 
to a wide range of applications such as the medical and industrial use of high 
intensity ultrasound in lithotripsy, thermotherapy, ultrasound cleaning, etc. 
There have been quite a few works in this aspect, more from engineering 
viewpoint; see [H da ESI la ED]. Rigorous study of the small global so¬ 
lutions and their decay estimates has been conducted in US]. A thorough 
study of the linearized models is a good starting point for better understand¬ 
ing the nonlinear models. Actually, the work m has shown that, even in 
the linear case, rich dynamics appear. 

Here we investigate the linear MGT system when memory effects get 
involved. For reasons that will be clear soon, we start with introducing 
wave equation. 

It is well known that wave system conserves mechanical energy. Gonsider 
the wave equation 

utt — Au = 0. 

Multiplying it by ut, which is called a “multiplier”, we have 

+ IIVull") = 0 ^ E{t) ^ i(||»,|p + IIVulp) = £(0). 

Here E[t) represents the mechanical energy, summation of kinetic and po¬ 
tential energy. We see the mechanical energy does not change when a wave 
evolves. On the other hand, we normally see mechanical energy decays in a 
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physical system due to certain damping mechanism. A different viewpoint 
is, in order to force the energy decay, we have to inject damping mechanisms 
into the system. There are different ways to implement this idea. One wajQ 
is to add viscous damping, also called frictional damping, into the wave 
equation. Taking into account the friction, we end up with equation 

utt - Au + ut = 0. 


The same multiplier ut gives energy estimate 

+ IIVulP) = -||u.|t < 0 ^ £(() < E{0). 

While the energy functional E{t) keeps the same form i(||ut|p + ||Vtt|p), the 
multiplier ut generates a pure negative term — ||ut|P on the right, causing 
damping in the system. A quick thought would be, if ||ut|p controls the full 
energy, say ||ut|p > CE{t), then we have ^E{t) < —CE{t) and Gronwall in¬ 
equality immediately gives the exponential decay of energy, E{t) < 

It is not the case here, however. 

Multiplier ut only gives us partial information: the mechanical energy 
is decreasing (at least when ut ^ 0). Without further investigation. It 
is not clear how fast it decreases, or whether it decreases exponentially. 
Fortunately, the energy here does decay exponentially. Indeed, it can be 
shown, by standard Lyapunov function method (essentially with the help of 
another multiplier u), that the energy decays exponentially [8]. 

What interests us here is a different kind of damping mechanism, namely 
the energy dissipation caused by viscoelasticity, which forces the appearance 
of memory term in the system. Viscoelasticity is the property of materials 
that exhibit both viscous and elastic characteristics when undergoing defor¬ 
mation. It usually appears in fluids with complex microstructure, such as 
polymers. One encounters viscoelastic materials in biological science, mate¬ 
rials sciences as well as in many industrial processes, e.g., in the chemical, 
food, and oil industries. The phenomena and mathematical models for such 
materials are more varied and complex than that of pure elastic materials 
or that of pure Newtonian fluids, see|27j. 

Our work is motivated by a series of papers on the wave equation (or its 
abstract version), with memory 


( 6 ) 


Utt — Au + / g{t — s)Au{s)d. 

Jo 


s = 0, 


where A is the Laplacian operator defined on certain functional space. The 
convolution term g(t — s)Au(s)ds represents memory: the integral itself 
suggests the nonlocality in time; the system at present moment “remember” 
information in the past history. Generally g(t) is assumed to be decreasing, 


^Adding structure damping —Aut to the wave equation is another way to obtain ex¬ 
ponential decay of the energy. These topics are of great importance in real applications, 
hence forming an active research area [T]. 
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suggesting the “nearer past” has stronger influence on the system than the 
“further past” does. This is a simplihed one dimensional wave equation with 
memory. Like all the physically meaningful partial differential equations, it 
is derived based on the balance laws of momentum and/or conservation of 
mass, together with a constitutive equation relating stress to strain (elastic¬ 
ity) or stress to rate of strain (fluids). For more information on the derivation 
of this 1-D system, we refer to m- 

As for wave equation, there are hundreds of papers studying either the 
damping effects caused by the memory terms or the interaction between dif¬ 
ferent damping mechanisms. See [SlEliaElEllTllIIllISllinilllllMllSSlIlH] to 
name a few. While a memory term does generate damping in a wave system, 
the nature of this mechanism seems to be subtle. In m, the authors found 
that a damped wave equation, which is already exponentially stable, loses 
this property after being added a memory that is not properly calibrated. 
This suggests, the more is not necessarily the better. It is important to un¬ 
derstand the potential benefits as well as drawbacks of adding memory to a 
system in applications. 

It is natural to ask questions such as, what kind of term creates damping 
effect in a system? is it possible that adding a damping term “hurt” the 
energy, in any sense? The answer hides in the energy equality in a wave- 
memory system. Roughly speaking, the memory term creates damping by 
weakening the total energy in the first place: it borrows energy to initialize 
the deformation and pays back in the long run. We exploit this point in 
another paper [ 22 ] . 

To our best knowledge, the present paper is the first work to investigate a 

temporally third order system that interacts with memory. As we will see in 

Section 2, the MGT equation is separated into different regimes based on a 

2 

critical parameter 7 = a — In the non-critical regime (7 > 0), which will 
be explained presently, a substantial portion of frictional damping exists. 
Hence it is not surprising that the energy in non-critical case intends to 
decay more often than that in critical case. 

A big part of the work is on how the memory term influences the MGT 
system in the non-critical case, namely for 7 7 ^ 0. This is not particularly 
challenging since the system alone contains strong damping mechanism. But 
this study brings up more interesting questions such as how different damp¬ 
ing mechanisms “cooperate” with each other. This is never a trivial issue, 
see miss]. Moreover, this step provides intuitions when it comes to the 
critical case 7 = 0 , where the system itself does not provide enough damp¬ 
ing. In the critical case, we show that “the right mixture” of memory can 
make the MGT system exponentially stabl^. Based on the results in wave 
system, this result is surprisingly neat, see Theorem 11.101 

The following notations are frequenctly used in the work. 


o 

It is currently an open question whether the first type of memory is sufficient to make 
the energy decay. 
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• A: a positive selfadjoint operator defined on a real Hilbert space 
H with domain 'D{A). We assume Aq > 0 exists such that ||u|| < 
Ao||^^/^u|| for all u ^ H. 

• ||u|| = ||u||jy: the norm on H. 

• -}{ = V{A^I^) X V{A^I'^) X H. 

• g o u = fg g(t — s)||u(t) — n(s)|p(is: an important quantity in the 
energy estimates. 

• g * u(t) = fg g(t — s)u(s)ds: the convolution integral. 

• G(t) = fg g(s)ds: the cumulative strength of viscoelasticity. 

• Fo(t) = + + the standard energy defined 

for MGT. 

• Fi(t) = Fg(t) — g' oA^^^u: the standard energy for MGT with mem¬ 
ory of type 1. 

• F 2 (t) = Fg{t) + g o A^^'^ut'. the standard energy for MGT with 
memory of type 2. 

• Fg{t) = Fg{t)+goA^^^w. the standard energy for MGT with memory 
of type 3, with w = Xu + ut- 

• Fg'"{t) = -|- ||rct|p + g o A^^'^w: energy for type 3 memory, 

in critical regime. 

• a,T,(?,b'- parameters in MGT equation. 

• ^ = a — a critical paremeter in defining concepts of non-critical 
and critical regimes. 

• Co,Cl, ..., etc: constants related to memory kernel g. 

• C, Co, Cl,...,etc: generic constants related to energy estimate. 

1.2. Main Results. We make the standing assumption 

Assumption 1.1. Assume 

• 1. g{-) G C^(0, oo) Pi C[0, oo), and g{t) > 0, g'{t) < 0 for t > 0. 

• 2. g''{t) > 0 for t > 0. 

• 3. Positive constants co,ci exist such that g' < —c^g. 

Remark 1.2. Among the above list, item 1 is the standing assumption on 
memory kernel. Item 2 says the kernel must be convex. This restriction can 
actually be removed; but to get energy bound, we need to impose condition 
on g(0). see [22]. Item 3 implies the exponential decay of the kernel. In 
|22j . we deal with problem that has memory kernel with general decay rate. 

Next we present the main results of energy decay in case of different mem¬ 
ory type. 

^Memory of type 1: 


We consider 


TUttt + (xutt + (?Au -|- bAut — / g{t — s)Au{s)ds = 0. 

Jo 
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We make the following assumption on the memory kernel. 

Assumption 1.3. Assume 

• 1 . 7 = 0 — ^ > 0 . 

• 2. G{+oo) < (?. 

Theorem 1.4. Under Assumptions [T^ and \l.tA the unique weak solution 
u exists on % and the standard energy Fi(t) decays exponentially, i.e., there 
exist positive constants C,uj > 0 such that 

Fi{t) < CFi{0)e-‘^\t> 0. 

^Memory of type 2: 


We consider 


TUttt + oiUtt + (?Au + bAut — / g{t — s)Aut{s)ds = 0 . 

Jo 


Assumption 1.5. Assume 

• 1 . 7 = 0 — ^> 0 . 

• 2. There exist positive numbers 6 and k £ f) that ^ < cq 
and 


G{+oo) < mm{ 


2{bk - c 2 ) 
2 + 9 


,b 



Theorem 1.6. Under Assumvtions [7771 and \1.5l the unique weak solution 
u exists on TL and the energy satisfies 


F2{t) < GF2{<d)e-^fit> 0 . 


^Memory of type 3: 


We consider MGT with a mixed memory in form of 

TUttt + oiUtt + (?Au + bAut — g * Aw = 0, 
where w = Xu{t) + ufit). 

2 

For the non-critical case, 7 = 0 — ^ > 0. Define the standard energy 
functional 

F^{t) = \\A^/‘^ut\\^ + \\A^/‘^u\\^ + \\utt\\^ +goA^/‘^w. 

Assumption 1.7. Assume 

• 1. 7 = 0 — ^> 0 . 

• 2. X £ (^, such that G{+oo) < 

Theorem 1.8. Accept the Assumvtions \1.1\ and [13 then the unique weak 
solution u exists on TL and the energy decays exponentially, i.e. 

F3{t)<GF3{0)e-‘^\t>0. 
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For the critical case where 7 = 0, the energy corresponding to the above 
system becomes 

A I |^l/2y;| |2 + I |2 + ^ o 

Assumption 1.9. Assume 
.F 7 = a-i^ = 0 . 

. S. A = f,G(+oo)<^. 

Theorem 1.10. Accept the Assumptions 1 1 . Il and \1.^ then the unique weak 
solution u exists on % and the energy F^'^{t) decays exponentially, namely 
we have the exponential decay o/||tct|| and ||A^/^u)||. 

Remark 1.11. The existence part of each proof shall be omitted; we focus 
on the study of decay rate. In fact the existence can be easily carried out by 
the standard Galerkin method. The main ingredient of Galerkin method is 
certainly the global a priori bounds of the solutions, which can be found in 
the stability results. 

Remark 1.12. In this work, we only consider exponential decay, which 
demands the exponential decay of the memory kernel, namely g' < —c^g. In 
paper m, based on the idea in PQ], we will deal with the more general case 
g' + H{g) < 0, where H{-) is a convex increasing function with II{0) = 0. 

The paper is organized as follows: In Section 2 we revisit the linearized 
Moore-Gibson-Thompson equation. The main result from |15] is stated and 
a slight different proof is presented. In Sections 3, 4 and 5, we give the 
proofs of the main theorems stated above. 

2. Revisit the MGT equation 

The MGT model is third order in time. The relaxation parameter r ac¬ 
counts for finite speed of acoustic waves, addressing the paradox of inhnite 
speed of propagation occurring in modeling acoustic waves [l5]. As a conse¬ 
quence of the presence of r, the original parabolic like model(r = 0 ) becomes 
hyperbolic(r > 0 ); the issues of well-posedness, stabilization and long time 
behaviors become more involved due to the presence of inhnitely many un¬ 
stable eigenvalues. With r > 0, roughly speaking when diffusion parameter 

6 = 0, ([5]) becomes ill-posed; on the other hand with b > 0 the system 
generates a strongly continuous semigroup. Moreover, the corresponding 
semigroup is exponentially stable on Ti when 7 = 0 : — (?Tb~^ > 0. When 

7 = 0 the semigroup is conservative. 

More specifically, consider the MGT equation ([5]), 

TUttt + autt + c^Au + bAut = 0, 

with initial conditions 


u(0) = Uo,Ut(0) = Ul,Utt{Q) = U2. 
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Here ^ is a self-adjoint positive operator on a Hilbert space H with a dense 
domain V^A) C H. The initial data 

{uo,ui,U 2 ) eH = X V{A^/^) X H. 

Here the parameters r>0,6>0,aG7^ represent the relaxation, the 
diffusivity and the friction. We state the following results from |15] : 

Theorem 2.1 (Kaltenbacher et.al. 2011 , [IS])- Let b > 0,a ^ TZ, then the 
system generates a strongly continuous semigroups on LL. 

2 

Moreover, introduce the parameter 7 = 0 — we have 
Theorem 2.2 (Kaltenbacher et.al. 2011, [IS]). Let 6 > 0,a > 0,c > 0. Let 

Ei{t) = b\\A^/'^{ut + + r\\utt + + yTikdP; 

E2{t) = a\\ut\\‘^ + 

and 

E{t) = Eiit) + E2{t). 

1. If "f > 0, then the semigroup generated in Theorem \2.1\ is exponentially 
stable on H. And there exist u > 0,C > 0 such that 

E{t) < Ce-^^E{<d)A > 0. 

2. If 'j = 0 , the energy Ei (t) remains constant. 

Remark 2.3. Theorems \2.1\ and \2.^ were proved in m- Here we give 
a slightly different proof of Theorem \2.2l In this process, the readers can 
familiarize themselves with the notations and the methodology that will be 
used in the later sections. 

Proof. 1. Non-critical case: 7 > 0 . 

We start with constructing energy functional and doing energy estimate. 

Multiplying (|5|) by uu, we have 

+ Roiit) = 0 
dt 

with 

Eoi{t) = [r||utt|p -h b\\A^/‘^ut\f + 2c^{Au,ut)], 

Roiit) = 2a\\utt\\‘^ - 
Multiply dSj) by ut to get 

-7-^02 (i) + Ro 2 it ) = 0 
dt 

with 

Eo 2 it) = [c^\\A^^^u\f + a\\ut\f + 2Tiuu,ut)], 

Ro2it) = 2b\\A^/\t\\^-2 t\\uu\\^. 
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Sinces 7 = a — ^ > 0, we can pick up a A; such that 


c , a 
— < k < —. 

b T 

Let -Eo(A) = -L'Oi(A) + kEQ 2 {i) be the natural energy functional, we have 

Eo{t) = b\\A^/‘^ut + + T\\utt + kut\f 

(7) 

+fcr(- - fc)||u,||^ + c^{k - 

T 0 


and 

( 8 ) 


^-® o ( a ) + Ro{t) — 0 , 


Ro{t) = 2T{--k)\\uuf + 2b{k-^)\\A^/^ut\\^. 

T b 


Remark 2.4. Note the coefficient of each term in Eq and Rq is strictly 

positive. This is possible only because of the existence of such a k with 
2 

\ < k < ^. This is critical to generating an effective damper in the system, 

as we will see soon. This is also why we call it “non-critical case” in this 

2 

work. For the “critical case”, namely when 7 = 0 — ^ = 0, no such a k 

exists. In other words, the k in the above calculations is forced to satisfied 
2 

k = Y So Ro{t) = 0 and we have no damper in the energy equation, 

hence no energy decay can be obtained this way. 


Remark 2.5. The natural energy functional comes directly from the 
energy equality and normally looks messy due to complicated expressions or 
additional parameters. So we get it dressed up and introduce the equivalent 

standard energy functional. 

We denote the standard energy functional by 

(9) Ffft) = WuttW^ + + WA^'M?. 

We have L^o(A) ^ Ffft), namely they are equivalent: there exists Ci > 
0 , C 2 > 0 such that 

CiEoit) < Fo{t) < C2Eoit). 

Indeed, recall ||Mt|| < this can be easily proved through the 

following lemma. 

Lemma 2.6. Let Cq > 0. Then 

||/ + 5|P + Co|blP~||/|P + lblP. 

Proof. Note 

11/ + 5 IP + C'olls'lP = ^ ll/lP + + (1 + -^)lblP 
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> Cidl/ll" + ll9lP), Cl = min{(l - 

The other direction \\f + gW"^ + Coils'll < <^2(11/11^ + 11(711^) is trivial. □ 
Similarly, for E{t) in Theorem 12.21 we can easily show 

E{t) ~ Tb(t). 

This simply means if we can show the exponential decay of Tb(t), we get 
the same decay for E{t). 

Generally, we want to prove the standard energy decays in certain way. 
Typically here we expect an exponential decay of Eq, namely, Fo{t) < 
Coe~^* for some positive constants Co,C. One way of doing this is to 
construct Lyapunov function L(t) based on the equation system, which is 
equivalent to energy To(t); then one expects the Lyapunov function to show 
up in form of Gronwall inequality, which guarantee proper decay. 

In the energy estimate, Gronwall inequalities take different forms. The 
most familiar one is in differential form, for non-negative function L{t) : 
[0, oo) —)■ R^, 

j^L{t) < -CL{t) => L{t) < L(0)e-^*. 

So if here we have 

< -CEoit), 

or equivalently 

< -cm, 

then we get the exponential decay. 

In our case, it is natural to pick L{t) = EQ{t). The problem is, from ([8|), 
we only have 

< -c{\\utt\\‘^ -k ^ -cEo{t). 

So the Gronwall inequality in differential form does not apply. We need it 
in more general form. The following lemma is taken from [21 p.l345]. 

Lemma 2.7. Let E{t) be a non-negative decreasing function on [0, oo) and 
if Cl > 0, To > 0 exists such that 

/ OO 

E{s)ds < CiE{t),t > To, 

then 

F{t) < T(0)exp(l - 0^4^), Vt > 0. 

Li -k -to 

From the energy identity ([8]), we know Eo^t) is non-increasing. Moreover, 
integrate over [t,T], we have 

Eo{T) + C3j^^\\uu\\^ + m/m?ds < Eo{t), 


( 10 ) 
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for arbitrary 0 <t <T and C 3 = min{2r(^ — k), 2b{k — ^)}. In particular 
we have 

WuttW"^+ \\A^^‘^ut\\‘^ds < ^Eo{t). 

Now if the following inequality holds for some constant 6*4 > 0 independent 
of t,r, 

^ \\A^/^ufds<CiEo{t), 

we would arrive at 



s < C 2 



s < CEq 


it) 


for arbitrary 0 < t < T < 00 . Hence take T = 00 , we would have the 
exponential decay of energy by Lemma 12.71 
To show 

\\A^/^u\\‘^ds < C^Eoit), 

multiply (l5|) by u, 



+ c^WA^/^uW^ = a\\ut\\^ + ^[^\\ut\f - T{utt,u) - a{ut,u)], 

which gives 

^\\A^^‘^u\\‘^\J+c^ \\A^^‘^u\\‘^ds = a \\ut\\‘^ds+[^\\ut\\‘^-T{utt,u)-a{ut,u)]\f 

^ ^\\A^^‘^u{T)\\‘^ + [ \\A^t‘^u\\‘^ds 

^ Jt 

< ^\\A^^'^u{t)f + a \\ut\fds + [^WutW^ -t{uu,u) - a{ut,u)]\J 

<a WutW^ds + CFoit) + CFo(r). 

Since ||ut|| < Ao||^^/^ut||, in view of (1101) . the equivalence of Tb(t) and Eo{t) 
and the non-increasing property of Eo{t), we have 

[ \\A^^‘^u\\‘^ds <C4Eo{t). 


So we get the exponential decay for £'o(t), hence for To(t) and E(t), since 
they are all equivalent. The proof of first part of Theorem 12.21 with 7 > 0 is 
completed. 

2. Critical case: 7 = 0 . 
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In the critical case, 7 = a — ^ = 0, all the calculations valid with 
k = (?/b = ce/r, but the expressions shrink because of cancellations. We 
have 

~ b\\J^^'^ut + —+ T||utt + 

= b\\A^/'^{ut + ku)\\^ + r||(ut + ku)t\\^, 

and 

j^E^it) + R^(t) = 0, 

B^it) = 2t(- - k)\\utt\\^ + 2b(k - d)||^i/2„,||2 ^ 0 . 
r 0 

Here E^[t) represents the energy functional in the critical case. R'^it) = 0 
here suggests the nonexistence of damping. That is, we have 

It is not difficult to see that EQ^{t) is exactly Ei{t) when 7 = 0 . So the 
energy is conserved. This completes the proof of Theorem 12.21 □ 

Remark 2.8. This is exactly the intuition provided by m- the MGT equa¬ 
tion IB can be written as 

7c^ 

TZtt + bAz + ■yzt = 

b 

2 

where z{t) = ut^u. Since r > 0,6 > 0, it is easy to see that, in the new 

variable z, becomes wave equation when 7 = 0 , hence no decay can be 

expected. This is exactly how MGT is connected to wave equation through 

2 

the parameter 7 = a — On the other hand, when 7 > 0, the left side 
of ^ becomes wave equation with frictional damping which intends to force 
exponential decay. 

Remark 2.9. From the spectrum viewpoint, this fact can he seen as well. 
Taking the prototype of ^ as an example where A = —A. After Fourier 
transform, letting u = F[u], we have a third order ODE in time 

TUttt + autt - b^^ut - c^i‘^u = 0 . 

It is not difficult to verify, the spectrum (roots of its characteristic equation 
rr^+ar^ ——= Oj are strictly negative as long as 'y = a — c^rlb > 0; 
while one of the roots becomes zero when 7 = 0 . 

From next section, we start to study MGT with memory terms. 
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3. Memory of type 1: Proof of Theorem 11.41 


We consider the equation ([5|) with memory term of the Hrst kind, so 
w{t) = u{t). This leads to 



( 11 ) 


Proof of Theorem 11.41 

Again we use multipliers and work through the energy estimate in several 
steps, where we define energy functional and show the decay of the solution. 
When multiplying (llip by uu and ut , only the memory terms need additional 
calculations because the other terms are the same as in last section. 

Note that here for the memory term, we have 
a.) inner product with uu- 



with 


(12) Eum — —g 



(13) 


= -g” o + g'{t)\\A^A u\\^ 


b.) inner product with ut- 



with 


(14) 


Ei2m = goA^/‘^u - G{t)\\A^Au\\^ 


(15) Ri 2 m = g' o A^^'^u - g{t)\\A^^'^u\\'^ 


Now we are ready for the energy estimate. 


Step 1. Multiply by uu, we have 


-;-£^ii(i) + Riiit) — 0, 

at 


with 


En{t) — Eoi{t) + Eiim{t), 

Ru{t) = Roi{t) + Rlm{t), 


where Enm, Riim are defined in (fT^ and (fT^ . Also recall 
Eoi{t) = [r||u4t|p + 2c^{Au,ut)] 


and 


Roi{t) = 2a\\uu\\'^ - 2c^||.4^/^Ui|p 
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as defined in the previous section. 

Step 2. Multiply by ut, we have 

-^Ei 2 {t) + Ruit) = 0 

with 

El2{t) = Eo2{t) + Ei2m{t), 

Rl2{t) = Ro2{t) + Rl2m{t)- 
Here Ei 2 m, Ri 2 m are defined in (fH)) and IfTS]) and 

Eo2{t) = [c^||M^/^u|p + Q;||uf|p + 2 r(uti,Ut)], 

Ro2{t) = 26||M^/^ut|p - 2T\\uu\f. 

2 2 

Step 3. Since 7 = 0 — ^ > 0, we can pick a k € (^, and define the 
natural energy Ei{t) = Eii{t) + kEi 2 {t). We have 

^Eiit) + Riit) = 0 
at 

with 

Ei{t) = £'o(t) + Eiim + kEi2m, 

Rl{t) = i?o(0 + -Rllm + kRi2m- 
Here Eq,Rq are from ([7|) and ([ 8 ]). 

Step 4. Now denote the standard energy 

(16) Fi{t) = lluttiP + - g' o 

We claim that 

Lemma 3.1. Ei{t) ~ Ei{t) 

Proof. From ©-(jS]), (fT^ - (fT5]) . we have 

Ei{t) = b\\A^^'^ut\\^ + 2(?{Au,ut) + (?k\\A^^‘^u\\^ 

+T\\utt + fcutlp + /cr(- - k)\\ut\\^ 

T 

— {g' — kg) o A^^'^u + {g — kG{t))\\A^^‘^u\\'^ 

+2 f g{t — s){A^^‘^u{t) — A^^'^u{s),A^^'^ut{t))ds 

Jo 

-2G{t){Au{t),ut{t)) 
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= + 2(c^ — G{t)){Au, ut) + k{c^ — G{t))\\A}^‘^u\\^ 

+T||itit + kutW^ + kT{- - k)\\ut\\^ 

T 

— {g' — kg) o 

+2 f g{t — s){A^^‘^u{t) — A^^'^u{s),A^^'^ut{t))ds. 

Jo 

= {b-^^^)\\A^/^ut\\‘^ 

+ -—+ 2k{Au,ut) + k‘^\\A^^‘^u\\'^ 

+t\\uu + kutW'^ + ^''■(- - k)\\ut\\'^ 

T 

— {g' — kg) o AJI'^u + 

+2 f g{t — s){A^^'^u{t) — A^^‘^u{s),A^^'^ut{t))ds 

Jo 

= ‘^- °(‘> ||^l/2„, + + (i, _ ^)||^'/2„,||2 

k k 

+T\\utt + kutW"^ + kT{— - k)\\ut\\‘^ 

T 

(17) -g'oA^/‘^u + g\\A^/‘^u\\^ 

+^\\A^/\^\\^ + kgoA^I\ 
k 

+2 f g{t — s){A^^'^u{t) — A^^'^u{s),A^^'^ut{t))ds. 

Jo 

Since 

2 | [ g{t—s){A^^'^u{t)—A^^'^u{s),A^^‘^Ut{t))ds\ < kgoA^^‘^u+^^\\A^^'^Ut\\'^, 

Jo k 

the summation of the last three terms is nonnegative. So 

El(t) > + kA^'^ul I" + (6 - ^)I112 

k k 

+T\\utt + kut\\^ + kT{^ - k)\\ut\\^ 

-g' o A^/‘^u + g\\A^/‘^u\\^. 

Recall g > Q,g' < 0,G{t) = g(s)ds < g{s)ds = G{+oo) < c^, by 
Lemma 12.61 we have 


Ei{t) > CFi(t). 
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On the other hand, it is easy to see 

= gw 11^1/2^^ + fe^V2^||2 + (5 _ i||^l/2^,||2 

k k 

+T\\utt + kutW^ + kr{- - k)\\ut\\‘^ 

T 

-g' o + g\\A^/‘^u\\^ 

+ M||yll/2„,||2 + tjo^'/2„ 

k 

+2 f g{t — s){A^^'^u{t) — A^^‘^u{s),A^^'^ut{t))ds 

Jo 

< C\\A^/\t + kA^^^u\f + C\\A^/^utf 

+^2*111111 + + g||^^t|P 

-g' oA^/^u + C\\A^/\\f 
+C\\A^^'^ut\\^ + Cg o A^/'^u 

< CFi{t) 

because of the boundedness of g{t),G{t) and the Assumption [T31 Hence the 
equivalence of ii'i(t) and -Fi(t) is established. □ 

Step 5. Recall 

Ri{t) = 2r(- - k)\\utt\\^ + 2b{k - ^)\\A^/‘^ut\\‘^ 

T 0 

(18) +{g'' — kg') o AJ^’^u + {kg — > 0. 

Since the right hand side does not contain a strictly negative multiple of 
we are one term away from the Gronwall inequality in derivative 
form. So we follow the idea in section 2, seeking the integral form. 

From the above energy equality we have 

(19) 

g J IkttiP + + ia" - kg') O A^^'^U + {kg - g')\\A^/'^u\\‘^ds 

+Ei{T) < Ei{t) 

for any 0 <t <T < oo. Note g" > 0, g' < 0, g > 0, in particular we have 

( 20 ) < C'£’i(t). 

The Gronwall inequality in integral form Ei{s)ds < CEi{t) would be 
established, if only we can show 

[ \\A^/^ufds<CEi{t). 
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Step 6. To achieve the goal, we apply multiplier u to (HD 


d h 
dt 2 




= - 'r{utt,u) -a{ut,u)] 

+G{t)\\A^^‘^u\\^ — f g{t — s){A^^'^{u{t) — u{s)),A^^‘^u{t))ds. 

Jo 

which gives 


(c^ — G{t) — eG{t))\\A^^‘^u\\^ds 


< a 


i: 


As + [-||ut|r - r{utt,u) - a{ut,u) 


\I 




-ll'lf. 


Since cP > G(+oo) > G{t), we can choose e so small that 


(c^ - G{t) - eG(t)) > > 0 


Since 


/ Wutlfds^Xlf WA^^'^utW'^ds < GEi{t), 

Jt Jt 

\\ut\?\J < + \\A^/^ut{T)\\^ < G{E^{t) + Ei{T)) < 

\{utt,u)\\ < \utt\\^+ < |nt|p < GEi{t), 

\{ut,u)\\ < GEi{t), 
we have ^ 

<GEi{t) + G goA^^‘^uds<GEi{t). 

The last inequality is because of (f20l) and the assumption g' < —cog, or 
equivalently g < ——g'h 

CQ 

T 

Step 7. Now combining all the calculations above, we have Fi[s)ds < 
GEi{t), hence 

[ Ei{s)ds<Gi [ Ei{s)ds <GEi{t). 

Jt Jt 

Applying Lemma 12.71 we have completed the proof of Theorem 11.41 □. 


Remark 3.2. Theorem \1.4\ is for non-critical case. For the critical case, 
the problem is open. Following all the calculations above, we arrive at 
but without terms in Ri like ||utt|p, as damper. In other words, 

only memory terms contribute to the damping. That creates challenging 
problems, which might require new multipliers or new methods. 
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4. Memory of type 2: Proof of Theorem 11.61 

We consider the equation ([5]) with memory term of the second kind, so 
w{t) = ut{t). This leads to 

( 21 ) TUttt + O-Utt + (?+ bAut — / g{t — s)Aut{s)ds = {) 

Jo 

Proof of Theorem 11.61 Next we work through the energy estimate. 

Step 1. Multiply by utt, we have 

d 


dt 


E2l{t) + R2i{t) — 0 


with 


and 


E 2 i{t) = [rlluitlp + + 2c^{Au,ut)] 

+goA^Ay _ G{t)\\A^Ay\\‘^ 


R 2 i{t) = 2a\\utt\\‘^ - 2 c^ 11 . 4 ^/^ 114 IP 
-g' o A^Ay + g{t)\\A^Ay\\ 2 _ 

Note v{t) = ut{t) here. 

Step 2. Multiply by ut, we have 


with 


and 


— E 22 {t) + R 22 {t) = 0 
dt 

R 22 {t) = 26||.4^/^n4|p - 2 r||n 4 i|p - 2G(t)||.4^/^u| 


+2 


[ g{t - s){A{v{t) - v{s)),ut{t))ds 
Jo 


E22{t) = [c^||.4^/^u|p + a||ii4|p + 2r(nt4,nt)]. 

2 2 
Step 3. Since 7 = 0 — ^ > 0, we can pick a k such that \ < k < ^. 

Let E 2 {t) = E 2 i{t) + kE 22 {t) be the natural energy, we have 

E 2 {t) = 6||.4^/^iii|p + 2(?{Au,ut) + c^A;||.4^/^m|P 


a 


+r||?iti + fciiill +A;r(- k)\\ut\\ 

T 

+goAA\-G{t)\\A^Ay\\\ 

= {b — G(t))||.4^/^ut|p + 2 c^(^tt, Ut) + c^A:||^^/^it| 
+T||n44 + feutlp + kT{— - A;)||nt|p 


and 

( 22 ) 


dt 


+g o AJAx 


E 2 {t) + R 2 {t) — 0 , 






MGT EQUATION WITH MEMORY 


19 


R 2 {t) = 2t(- - k)\\utt\\^ + 2 b{k - ^)\\A^/‘^ut\\^ 

T 0 

-g' o + g{t)\\A^/\\\^ 

—2kG{t)\\A^^'^v\\'^ + 2k f g{t — s){A{v{t) — v{s)),ut{t))ds. 

Jo 

Step 4. Now denote 

F2it) = \\uu\\‘^ + \\A^^‘^ut\\‘^ + WA^/^uW"^ + goA^/'^v, 
we claim that 

Lemma 4.1. E 2 {t) ~ -^2(0 

Proof. Recall G{t) = g{s)ds and G(+oo) < ci < b — 

E2{t) = (6 — G(t))||.4.^/^it4|p + 2 c^(.4,m, tit) + c^A:||.4.^/^u|p 

+ T\\utt + kutW^ + kri-- k)\\ut\\^ 

T 

+ go AJI‘^v. 

1 ^—■ 2 

By Assumption 11.51 Gft) < G(+oo) < 6 — we have 

{h — G{t))\\A^^‘^ut\\^ + 2(?{Au, ut) + (?k\\A^^‘^u\^ > c(||A^/^iit||^ + ||A^'^^u|p) 
for some positive constant c. So we have 

E2{t) > cF2{t). 

The other direction is trivial. So the equivalence of E 2 {t) and F 2 {t) is 
established. □ 

Step 5. Recall 

R 2 {t) = 2t(- - k)\\utt\\^ + 2 b{k - ^)\\A^/'^ut\\^ 

T 0 

oAJI’^v + g{t)\\A^/‘^v\\‘^ 

—2A:G(t)||A^'^^i;|P + 2/c f g{t — s){A{v{t) — v{s)),ut{t))ds 

Jo 

= 2r{^ - k)\\uu\\^ + 2k{b -j- Gm\A^/^utf 
-g' o AJI'^v + g{t)\\A^^‘^v\\‘^ 

+2k I g{t - s){A{v{t) - v{s)),ut{t))ds 

Jo 

Under the assumption 

g' < -cog => -g' > cog 

2 k f g{t - s){A{v{t) - v{s)),ut{t))ds 
Jo 


Now 
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<2k\^a{t)\\A'l\f + ^goA'l\\ 

= kea[t)\\A^i'‘v\f + 

u 

Under Assumption 11.51 a straightforward calculation shows 
R2{t) > c{\\utt\\^ + + goA^/‘^v) 

and 

(lluttll^ + + g o A^/‘^v)ds < cE2{t). 

The Gronwall inequality in integral form 

E{s)ds < cE{t) 

would be established, if we can show 

\\A^/‘^ufds < cE{t). 

And this is achieved by same idea as the previous section. We omit the 
calculation. So the Theorem [UHl is completed. □. 

5. Memory of type 3: Proof of Theorems 11.81 and 11.101 

Now we are concerned with MGT with mixed memory. The equation 
reads 





(23) TUtu + autt + c^Au + bAut 
where w{s) = \u{s) + ut{s). 


f 


g{t — s)Aw{s)ds = 0, 


Proof of Theorem 11.81 


Let E 3 ,R 3 denote the energy and damper in this case. Recall how A is 
defined. We pick k = X and combine the calculations in the first two types 
of memory to get 

-J-E 3 + Rg = 0. 

dt 

Here we have: 


( 24 ) 


+T\\utt + fcutlp + kT{- - A:)||ut|p 

T 

+g o 

Rg = 2(a - kT)\\utt\\^ + 2{bk - 
+g \p - g' o A^^'^w. 
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Let 

(25) Fs{t) = +goA^/‘^w. 

Now we prove the Theorem II.81 namely, F^{t) decays exponentially. 


Proof. Firstly, it is easy to see Fl 3 (t) ~ Fsit)- Also since E'^{t) = —R 3 {t) < 0, 
we have F' 3 (t) non-increasing. Moreover, 

C IkttiP + WA^^'^utW'^ + g\\A^^‘^w\\‘^ + g o A^^^w{s)ds 

< R^{s)ds < Esit). 

In particular, 

Il'WitlP + \\A'-^‘^ut\f + g o A^^'^w{s)ds < CE‘i{t). 

Now if we could show 

(26) J < cE3{t), 



we would complete the proof of Theorem 11.81 by Lemma 12.71 following the 
same idea as before. 

To show (f26|) . multiply (f23]l by u, we have 


0 = T—{utt,u)-T{utt,ut) 


+a—{ut,u) - a||Mt|p 
dt 




b d 
2 dt 


\\A^/\\ 


— / g{t — s){Aw{s),u{t))ds 
Jo 


dt 


2 dt 


d 


+a^{ut,u) - a||Mt|p 

+ c 2 ||. 4 '/ 2„||2 + ^|||^ 1 / 2„||2 

+ [ g{t — s){A{w{t) — w{s)),u{t))ds 
Jo 


-XamA^l^uf - id|G(t)||A/2„||2] + l9(t)||^i/2„||2. 
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That leads to 


(c^ - XG{t) + ^g{t))\\A^/^u\f 
d T d 2 

+e| \A^/‘^u\\^ + Ceg o A^/'^w 


+ 


dt 2 


Hence, 


i: 


(c^ — XG{t) + ^g{t) — e)\\A^^‘^u\\^ds 


< GE 3 {t) + Cj (WutW"^ + g o A^^'^w{s))ds 

< GEsit). 

2 

Recall G(+oo) < choose e sufficiently small and t sufficiently large so 
that —XG{t) — g{t)—e is strictly positive, we get ([26]) . hence the exponential 
decay of E^^t). The same for E 3 {t) since they are equivalent. □ 

Proof of Theorem 11.101 

We are considering 

TUttt + c^utt + c^Au + bAut — g * Aw = 0, 

with 7 = 0, A = ^. We shall prove + ||rct|p + go A^l'^w 

decays exponentially. 


Proof. Impose the conditions k = ^ into 
for the critical case 


, we get the energy estimate 


with 


|Er + flr = o 


Ef = 




+T||r(;ilp + go A^^'^w. 

— g' o A^l'^w. 

2 

It is easy to see E'^ is equivalent to since G(+cx)) < 


Bf{ = 
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Because 7 = 0, as suggested in Remark 12.81 the equation can be written 
into 

(27) TWtt + bAw — g * Aw = 0 . 

This is exactly the problem we dealt with in Lasiecka et.al.[ 2 Tj. From 
analysis of the second order equations with memory, it is known that 
system {w,wt) € D{A^^'^) x H decays exponentially. Adopting result, 
have shown Theorem 11.101 
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